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The running QCD coupling in the pre- asymptotic region 
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We study deviations from the perturbative asymptotic behaviour in the running QCD coupling by analysing 
non-perturbative measurements of a s (p) at low momenta (p m 2 GeV) as obtained from the lattice three-gluon 
vertex. Our exploratory study provides some evidence for power corrections to the perturbative running propor- 
tional to 1/p 2 . 
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1. INTRODUCTION 

The standard procedure to parametrise non- 
perturbative QCD effects in terms of power cor- 
rections to perturbative results is based on the 
Operator Product Expansion (OPE). In this 
framework, the powers involved in the expansion 
are uniquely fixed by the symmetries and the di- 
mension of the relevant operator product. The 
above picture has recently been challenged 0- 
9, when it was pointed out that power correc- 
tions which are not a priori expected from OPE 
may in fact appear in physical observables. Such 
terms may arise from (UV-subleading) power cor- 
rections to a s (p), corresponding to non- analytical 
contributions to the /3-function. Clearly, the ex- 
istence of OPE-independent power corrections, if 
demonstrated, would have a major impact on our 
understanding of non-perturbative QCD effects 
and would affect QCD predictions for several pro- 
cesses. For example, ^ contributions may be rel- 
evant for the analysis of r decays [f||l]] . 

It would be highly desirable to develop a the- 
oretical framework where the occurrence of these 
effects is demonstrated and estimates are ob- 
tained from first principles QCD calculations. 
The results in |||[ can be considered as a first 
step in this direction: some evidence for an unex- 
pected ^ contribution to the gluon condensate 
was found by means of lattice calculations. 



The aim of the present work is to test a method 
to detect the presence of power corrections in the 
running QCD coupling. Non-perturbative lattice 
estimates of the coupling at low momenta are 
compared with perturbative formulae. The final 
goal is to investigate the conjecture that OPE- 
independent power corrections to physical observ- 
ables are linked to power terms in the running 
coupling. Although at this stage our work is ex- 
ploratory in nature and further simulations will 
be required to obtain a conclusive answer, our 
analysis provides some preliminary evidence for 
power corrections to a s (p) for a particular defini- 
tion of the coupling. 

The paper is organised as follows: in Section |^ 
we briefly review some theoretical arguments in 
support of power corrections to a s (p), illustrat- 
ing the special role that may be played by ^j- 
terms. In Section || we analyse the lattice data 
and present some preliminary evidence for power 
corrections. Finally, in Section || we draw our 
conclusions. 

2. WHY POWER CORRECTIONS? 

Power corrections to a s (p) can be shown to 
arise naturally in many physical schemes 
Such corrections cannot be excluded a priori in 
any renormalisation scheme. Clearly, the non- 
perturbative nature of such effects makes it very 
hard to assess their dependence on the renormal- 
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isation scheme. A term of order A 2 jp 2 is a strong 
candidate for a power correction to a s (p) . To see 
why, consider the interaction of two heavy quarks 
in the static limit and in the one-gluon-exchange 
approximation (for a more detailed discussion see 
||). The static potential V(r) can be written as 

V{r) oc a s Jd 3 p^^. (1) 

If one inserts in the above formula a running cou- 
pling of the form a s (p 2 ) » A 2 /p 2 , this results in 
a linearly confining potential. Similarly, consider 
the "force" definition of the running coupling: 

where again V(r) represents the static interquark 
potential. A linear confinement term in V(r) gen- 
erates a 1/Q 2 contribution to the coupling, whose 
order of magnitude is given by the string tension. 
This can be interpreted as a clue for the existence 
of a pp contribution, providing an estimate for 
its expected order of magnitude, at least in one 
(physically sound) scheme. Finally, power correc- 
tions to a s (p) also emerge if one assumes that the 
singularities appearing in the pcrturbative formu- 
lae for the running coupling are "removed" by 
non-perturbative effects ||. 

3. ANALYSIS AND RESULTS 

We shall compare non-perturbative lattice data 
for a s (p) with simple models where a power cor- 
rection term is added to the perturbative formula 
at a given order. The first problem is the possible 
interplay between power corrections and our igno- 
rance about higher orders of perturbation theory. 
In particular, for the scheme that we will con- 
sider, the three-loop coefficient of the /3-function 
is not known. Knowledge of such a coefficient 
would allow a more reliable comparison of our es- 
timates for the A parameter in our scheme with 
lattice determinations of A in a different scheme, 
for which the three- loop result is available JTo| . 
In fact, although matching the A parameter be- 
tween different schemes only requires a one-loop 
computation, the reliability of such a comparison 
rests on the assumption that the value of A in 



each scheme is fairly stable with respect to the 
inclusion of higher order terms in the definition 
of A. In practice, when working at two- or three- 
loop order, the value of A is still quite sensitive 
to the order of the calculation. Even within such 
limitations, we will argue that it is possible to es- 
timate the impact of three-loop effects and that a 
description with power corrections seems relevant 
even at that order. 

3.1. Choice of the coupling 

We need to measure a s (p) at low momenta 
(where power-like terms may be sizeable) and in 
a relatively wide momentum range. For this pur- 
pose, the best method is one where a s (p) can be 
measured for several momentum values from a 
single Monte Carlo data set. One suitable method 
is the determination of the coupling from the 
renormalised lattice three-gluon vertex function 
0,0. By varying the renormalisation scale p, 
one can determine a s (p) for different momenta 
from a single simulation. Obviously the renor- 
malisation scale must be chosen in a range such 
that finite volume effects and discretisation errors 
are both under control. The numerical results for 
a s (p) used in this work were obtained by applying 
such a method on a sample of 150 Monte Carlo 
configurations on a 16 4 lattice at f3 = 6.0. The 
calculation was performed in the Landau gauge. 
For full details of the method we refer the reader 
to Ref. where such results were first pre- 

sented. In order to detect violations of rotational 
invariance, different combinations of lattice vec- 
tors have sometimes been used for a fixed value 
of p 2 , which accounts for the graphical "splitting" 
of some data points. 

3.2. Two-loop analysis 

At the two-loop level, we consider the following 
formula: 

r x = 1 _ h log(log( P 2 /A| ; )) 

s[P> 6„log(p 2 /A|) b (b \og(p 2 /A 2 2l )) 2 
A 2 

+ c 2l ±f (3) 
p z 

By fitting the data to (|J) we obtain two es- 
timates for (A 2 /,c 2 ;), namely (0.84(1),0.31(3)) 
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and (0.73(1), 0.99(7)), with comparable values for 
Xdof — The momentum range for the fit cor- 
responds to p ~ 1.8 — 3 GeV. We take the first set 
of values as our best estimate of the parameters 
as the corresponding value of A21 is close to what 
is obtained from a "pure" two-loop fit, i.e. A21 is 
stable with respect to the introduction of power 
corrections. This choice will be supported also 
by independent considerations at the three-loop 
level. In summary, a two-loop description with 
power corrections based on (||) fits well the data 
in a consistent momentum range. Our best fit is 
shown in Figure 1. 




Figure 1. The best fit to (g). 



3.3. Three-loop analysis 

A major obstacle for a three-loop analysis is 
our ignorance of the first non-universal coefficient 
&2 of the perturbative /3-function. In order to 
gain insight, we perform a two-parameter fit to 
the "pure" three-loop formula, taking and the 
unknown coefficient 62 as the fitting parameters. 
We call b^^ the fit estimate for 62- We obtain 
A 3/ = 0.72(1), b\ fi = 1.3(1), with xlof « 1-8 
(see dashed curve in Fig. 2). The momentum 
range where we obtain the best description of the 



data is p ~ 2 — 3 GeV. Our result for A 3 ; provides 
(via perturbative matching) an estimate for Aj^g 
which is in very good agreement with the estimate 
in |ict| , which was obtained from the computa- 
tion of the A parameter in a completely different 
scheme. Although our estimate depends on the 
extra parameter , the agreement between the 
two results is remarkable. 

So far, the success of the "pure" three-loop fit 
suggests that the power term in the two-loop for- 
mula merely provides an effective description of 
three-loop effects. However, it turns out that 
there is room for a power correction even at the 
three-loop level. To see this, we consider a three- 
loop formula with a power correction: 



1 _ h log(L) 
b L b {b Q L)i 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 



Figure 2. Fits to ([!]) (solid line) versus a pure 
three-loop fit (dashed line). 
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where L — log(p 2 /A 2 ; ) and is again to 
be determined from a fit. Fitting the data to 
§), we obtain A 3i = 0.72(f), b% ff = f.0(f) and 
c 3 = 0.41(2), with Xdof ~ 1-8, i n a momen- 
tum range 1.8 GeV < p < 3GeV (see Fig. 2). 
We note that the value for A3/ is fully consistent 
with the previous determination from the "pure" 
three-loop description. The value for b\^ is also 
reasonably stable with respect to the previous de- 
termination. By comparing results from fits to (^) 
and ([!]), it emerges that 



c 2 A 2 2l = 0.22(2) GeV 2 



c 3 A| 



0.21(2) GeV 2 . (5) 



In other words, the power terms providing the 
best fit to (||) and ([|) are numerically the same, 
so that there seems to be no interplay between 
the indetermination connected to the perturba- 
tive terms and the power correction term, within 
the precision of our data. We take this fact as an 
indication that a description in terms of power 
corrections is still relevant at the three-loop level. 
Notice that the numerical value of the power cor- 
rection is comparable to the standard estimate for 
the string tension. 

One could object that at the two-loop level 
we had chosen between two sets of values for 
(A2/,C2/), and that our choice is crucial for the 
validity of (Q). An a posteriori justification for 
our choice is obtained from the following test: 
we plot a few values for a s (p) as generated by 
the "pure" three-loop formula for A31 — 0.72 and 
62 = 1-0. Then, by fitting such points to the 
"pure" two-loop formula, one gets A21 ~ 0.84, i.e. 
the value for which (|^) holds. Again, the above 
test seems to confirm that perturbative and non- 
perturbative contributions do not mix in our for- 
mulae when upgrading from a two-loop to a three- 
loop description, thus suggesting that a genuine 
■p- correction is present in the data. 

4. CONCLUSIONS 

We have discussed an exploratory investigation 
of power corrections in the running QCD coupling 
a s (p) by comparing non-perturbative lattice re- 
sults with theoretical models. Some evidence was 
found for 1/p 2 corrections, whose size would be 
consistent with what is suggested by simple argu- 



ments from the static potential. 

Our results need to be further tested by the 
analysis of a larger data set and by a study of the 
dependence of the fit coefficient on the ultraviolet 
and infrared lattice cutoff. A very delicate issue 
is the assessment of the scheme dependendence 
of our results. In particular, we note that the 
definition of the coupling that we adopted is a 
priori gauge-dependent. This point will be the 
focus of our future work. 
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